Electrostatic interactions in critical solvents 
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Abstract 

The subtle interplay between critical phenomena and electrostatics is investigated by considering 
the effective force acting on two parallel walls confining a near-critical binary liquid mixture with 
added salt. The ion-solvent coupling can turn a non-critical repulsive electrostatic force into an 
attractive one upon approaching the critical point. However, the effective force is eventually 
dominated by the critical Casimir effect, the universal properties of which are not altered by the 
presence of salt. This observation allows a consistent interpretation of recent experimental data. 

PACS numbers: 05.70.Np, 05.70.Jk, 61.20.Qg 
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Introduction. Effective interactions among surfaces in contact with fluid media play a 
central role for a variety of topical fields in soft and condensed matter physics, cell biology, 
colloid and surface science, and nanotechnology. Since many relevant fluid media contain 
polar liquids such as water, their confining surfaces acquire an electric charge due to ion 
association or dissociation. As a result, electrostatic forces are expected to contribute signif- 
icantly to the interactions. In addition, a strong and highly temperature-sensitive solvent- 
mediated effective force arises upon approaching critical points. This critical Casimir force 
has recently been reported for a single colloidal particle close to a wall and immersed in 
a binary liquid mixture near its critical demixing point This force is expected to play 
also a role for the aggregation of colloidal suspensions [2|, yj. Motivated by recent exper- 
iments, here we investigate the interplay between electrostatic and critical Casimir forces, 
which turns out to be responsible for rather unexpected effects in binary liquid mixtures 
with added salt. Colloids dispersed in such a solvent have been reported Q to aggregate 
at a temperature difference from the critical demixing point which increases upon increas- 
ing the ionic strength, i.e., the screening of the electrostatic forces. This observation has 
been confirmed^experimentally also for a single colloid near a wall [4]. While it was orig- 



inally argued [2j that the aggregation could be completely explained in terms of a simple 
superposition of the critical Casimir and screened Coulomb forces (see, however, Ref. |3|), 
subsequent experimental results challenge this picture: An attractive colloid-wall interaction 
has been observed within a suitable temperature range even though both the electrostatic 
and the critical Casimir force are expected to be separately repulsive under these experi- 
mental conditions Q] . This points towards an important and yet unexplored aspect of the 
coupling between electrostatics and the critical fluctuations of the medium. Certain features 
of ion-solvent coupling near critical points were investigated in the past, such as the possi- 
bility of a micro-heterogeneous phase (sf] and the influence of criticality onto the solubility 
of ions 6|. However, the complementary point of view, i.e., the influence of ions onto the 
critical fluctuations of a solvent and therefore onto the critical Casimir effect has not yet 
been studied. We present a minimalistic but sufficiently enriched theory which explains the 
aforementioned unexpected experimental results. This contribution is intended to initiate 
a cross fertilization between research areas which so far focused separately on the critical 
Casimir effect in salt-free systems or on electrostatic interactions in fluctuation-free solvents. 
We expect that understanding — and thus being able to use on purpose — the coupling 
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between the critical Casimir effect and electrostatics provides a key to push forward the 
topical fields mentioned above. Our findings actually reach beyond the realm of soft matter 
because the counterintuitive behavior en route to universality discussed here is expected to 
be paradigmatic for many branches of physics in which generically critical and non-critical 
fields are coupled. 

Model. In a three-dimensional (d = 3) Cartesian coordinate system, we consider two 
parallel walls located at positions z = and z = L > 0, respectively, with the space in 
between being filled by a binary liquid mixture. In this solvent cations (+) and anions 
(— ) are dissolved. The solvent particles are all assumed to be of equal size with volume 
a 3 > whereas the ions are considered point-like. At the dimensionless position z := zfa, 
the number densities of the solvent components are given by <j)(z)a~ 3 and (1 — <j)(z))a~ 3 
with < <fi < 1, whereas the densities of the cations and anions are given by g + (z)a~ 3 and 
Q-(z)a~ 3 , respectively. The walls carry surface charge densities aoed~~ 2 at z = and cr^ea -2 
at z = L := L/a, where e is the elementary charge. The composition (f>(z) couples to surface 
fields ho at z = and Hl at z — L, where foo,z > (< 0) leads to a preferential adsorption of 
the solvent component with = 1 (=0). The equilibrium profiles 0, g + , and Q- minimize 
the approximate grand potential density functional fceTfi[0, g±], 
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with w so i(0) = 0(ln0-^) + (l-0) ln(l-0)+x(T)0(l-0) and u^(g ± ) = Q±Qng±-l-n±) 
as the grand potential bulk densities of the solvent and the ±-ions, respectively. Here k B T 
is the thermal energy, Aa 2 is the area of one wall, fi^hsT and fi±k B T are the chemical 
potentials of the solvent composition and the ±-ions, respectively, and = e 2 / (AirekBT) is 
the Bjerrum length for a uniform permittivity e; a ^-dependent permittivity j?J| corresponds 
to modified surface fields ho t L S|. The (temperature-dependent) Flory-Huggins parameter 
x(T) > describes the solvent-solvent interaction, which leads to a phase separation in 
the range x{T) > x(T c )] the gradient term oc <fi'{z) 2 accounts for the spatial variation 
of the solvent composition 9j. The ion-solvent interaction is given by the effective ion 



potential /csTV±(0) due to the solvent, with V±(0) = — ln(l — 0(1 — exp(— /±))) and where 
f±ksT is the difference of the bulk solvation free energies of a ±-ion in solvents with 0=1 
and = 0. This expression of V±(0) leads to a bulk phase diagram with a critical point 
which is shifted towards larger values of x{T c ) as the salt concentration increases V±(0) 
is an improvement of the standard approximation (j)f± (which it reduces to for /± <C 1) 



is 



which are, however, not 



because for f± > 1 the latter leads to multiple critical points 
observed experimentally. For f± — > oo the ion potentials V±(0) reduce to — ln(l — 0) which 
describes the entropy loss and thus free energy increase due to the insolubility of ions in 
the solvent component with = 1. The electric displacement D(z, [g±])ea~ 2 in Eq. (TTj) 



fulfills Gauss' law with fixed surface charges [10j: D'(z, [g±\) = Q+{z) — g_(z), D(0, [g±]) 



cr 0) D(L, [g±]) = —&L- Note that D(z,[g±\) is generated by the ±-ions and the surface 
charges cx^l, independent of 0. Within the present model, ions interact with the walls only 
electrostatically. 

In order to highlight the effect of the ion-solvent coupling we focus on an approximate 
grand potential functional for the solvent composition alone, which is obtained by expand- 
ing f2[0, g±] in Eq. (CQ) in terms of the order parameter tp := — 0& and the ion density 
differences A^ := g± — I retaining quadratic contributions as well as terms proportional 
to (f 3 and (p 4 . Here 0& and / = g±b denote the bulk solvent composition and the bulk ionic 
strength, respectively, corresponding to the chemical potentials \i§ and \i±. The minimiza- 
tion with respect to Ag± leads to linear, analytically solvable Euler-Lagrange equations for 
Ag± ccl (z, [if]), which are functionals of ip. Inserting these solutions into Eq. (TjQ) one obtains 
a Ginzburg-Landau-type functional 

m = j Az s xU(zMz) + m v{z? 



+ ^) 4 + ^'W 2 } (2) 
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with the temperature-like variable t(T) := 1/06 + 1/(1 — 4>b) ~~ 2x(7 1 ). Here we assume that 
the mixture is at its critical composition such that there is no 9? 3 -term. The electrostatic 
effects are contained in the coupling g := 2/0^ + 2/(1 — b ) 3 + 6/(7^ + 7i),7± := V±(4>i>), 
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as well as in an "external" field generated by the surface charges cto,l : 

{Z> ' 2(1 - exp(-2«L)) 

x [(o- + a L exp(-KL)) exp(-Kz) + (3) 

(cr L + do exp(-KL)) exp(-«(L - z))] 

with the Debye screening length k _1 = (87t£b/) _1 / 2 and A7 := 7 + — 7_. The "direct", 
i.e., solely ion-mediated, electrostatic interaction between the walls is given by W(L) := 
(4ii£b/ k)P{kL, (To, (Jl) where 

P{XML) := gM£±M±^Eizf) . (4) 

2 smh(x) 



The ion-solvent coupling affects the critical point only at order 0((A'~f) 2 ) 

Upon approaching the critical point the dimensionless bulk correlation length £ = £/a 



■\/ x(T) / (3t(T)) , which characterizes the exponential decay of the two-point correlation func- 
tion, diverges. Accordingly, on the scale £, U(z) is localized at the boundaries and therefore 
it merely modifies the surface fields /io,L- Consequently H turns into a standard y9 4 -theory, 
which describes the critical behavior of the Ising universality class Thus, within the 
present model, electrostatic interactions do not affect the universal critical behavior of the 
solvent. 

The effective wall- wall interaction is defined by 00 (L) := u(L)kBTa~ 2 with u(L) : = 
(H(L) — W(oo))/A, where H(L)/A is the minimum of Eq. (j2J). In general, for the critical 
contribution one has u(L) = q 9{L/^)/L d ~ 1 with a universal scaling function $(x), which 



depends only on the relative signs of /iq,l 12|, with $(x — > 0) = const and $(x — >• 00 



Cx d ~ x exp(— x), where C is a universal, boundary-condition-dependent constant 12]. 

For a sufficiently small bulk correlation length, i.e., if £ln£ C L, the term oc ip A in Eq. (J2J) 
can be neglected relative to the term oc (f 2 . The resulting quadratic functional can be readily 
minimized and leads to the approximate effective wall-wall interaction 

co(L) = - JLp(L/£, ho, h L ) + ^P(kL, * , a L ) 
X{T) k 

-A'y—-—(Qi(KL,K^)(h Q a L + h L a Q )+ (5) 



Q 2 (kL, Kt){ho<T + h L a L )) + C((A 7 ) 



2\ 



with the function, which is analytical for y > 0, 



n ( \ _J ( yexp(-kx/y) exp(-kx) . 

W , 2/J • y2 _ ! 1, ! _ exp(-2x/y) 1 - exp(-2x) , j ' 



In Eq. ((SJ) the term oc P(L/£, ho, hi) corresponds to the contribution of the surface fields 
to the wall-wall interaction in the absence of ion-solvent coupling (A7 = 0) whereas the term 
oc P(kL, o"o, (Tl) is the direct electrostatic wall-wall interaction. The term oc (h o~i + hicr ) 
is the interaction (oc A7) of the order parameter profile close to one wall (oc /&o,l) with the 
polarization of the diffuse ion layer due to the surface charge on the opposite wall (oc ctl,o)- 
This coupling between the order parameter and the ion density profiles is the central result 
of the present analysis and has important consequences (see below). The analogous term 
oc (h Q ao + hioi) is small. 

Discussion. In the following we discuss the experiments with colloids alluded to in the 
Introduction. The predictions of the present model for two walls can be readily translated 
into those for the wall-sphere and sphere-sphere geometry by means of the Derjaguin ap- 
proximation, which is applicable at separations much smaller than the sphere radii 10] . It 



turns out that assuming additivity of Casimir and Coulomb forces, i.e., independence of the 
order parameter from electrostatics, is in general insufficient to explain the experimental 
observations, whereas the present model, which includes ion-solvent coupling, leads to a 
consistent picture. 

First we consider symmetric boundary conditions, (ho, hi) = (—,—), for which the ion- 
solvent coupling is masked by the strong direct electrostatic repulsion. This situation has 
been investigated experimentally with a suspension of hydrophilic spherical colloids in a 
water-oil mixture Q as well as with a single hydrophilic colloidal sphere in a similar water- 
oil mixture near a hydrophilic glass wall {4]. In the presence of salt aggregation (2} or strong 
wall-sphere attraction jj] has been observed upon approaching the critical point of the binary 
mixture already several Kelvin away from the critical point. Within the present model, this 
setting is described by ho = hi < and a = ol < with the composition <p expressed as the 
mole fraction of the non- aqueous component. For a certain choice of parameters Fig. Ufa) 
displays the effective wall- wall interaction potential u(L). Since A7 ^ corresponds to 
/+ ^ /_, a negative ion-solvent coupling strength A7 < describes a salt the cations of 
which are slightly better soluble in oil than the anions, which is expected because the oils 
used in the experiments, 3-methylpyridine and 2,6-dimethylpyridine (2,6-lutidine), are Lewis 



bases 



13J. The parameters used in Fig. [T] correspond to a critical water-2,6-lutidine mixture 



(a = 0.34nm, £3 = 2.82) with lOmM salt (1/k — 7.73). Far away from T c the effective wall- 
wall potential oo(L) exhibits a repulsion due to the direct electrostatic wall- wall interaction. 
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Upon approaching T c , u(L) starts to develop an increasing attraction due to the critical 
Casimir effect. Since the change from repulsion to attraction occurs at « 1, the attraction 
sets in only very close to the critical point if the ionic strength is small, whereas this change 
occurs already considerably far away from the critical point if the ionic strength is large. 
Due to the strong direct electrostatic wall-wall interaction between hydrophilic walls, the 
ion-solvent coupling does not qualitatively influence the effective wall-wall potential, so that 
the assumption of additivity of critical Casimir and Coulomb forces [2, s| is justified. 

The situation is different for antisymmetric boundary conditions, (ho, hi) = (+,—), as 
studied experimentally in Ref. (4] using a single hydrophilic colloid near a hydrophobic glass 
plate. Repulsion is observed far away from as well as close to the critical point, whereas 
within an intermediate temperature range a strong attraction is found. The near-critical 
repulsion is readily understood in terms of the critical Casimir effect for antisymmetric 
boundary conditions and the repulsion far away from the critical point is of electrostatic 
origin. However, the attraction occurring in the intermediate temperature range cannot 
be explained within a picture without ion-solvent coupling. Figure DJb) shows oo(L) for a 
particular choice of non-symmetric surface fields ho > 0, hi, < and surface charge densities 
ct ,<Tl < 0, |o"o| <C |cr L | corresponding to a weakly charged hydrophobic wall. Far from 
( K £ < 0.28) and close to (k£ > 0.82) the critical point u(L) is repulsive because in Eq. ([5]) the 
terms oc P dominate. Upon increasing k£ beyond 0.28, i.e., en route towards T c , attraction 
occurs (see = 0.39), which for the chosen parameters is strongest around k£ = 0.46 and 
which weakens again closer to T c (see k£ = 0.59). This attraction is caused by the coupling 
(A7 < 0) between the order parameter profile near the hydrophobic wall (oc ho > 0) and the 
electrostatic potential due to the opposite, hydrophilic wall (oc oi < 0). For this effect to 
take place it is essential that the hydrophobic wall is sufficiently weakly charged ( | cr 1 <c \ctl\, 
see Ref. Li). 

n 

Ion-solvent coupling manifests itself in yet another experiment described in Ref. |4J, in 
which the surface preference of the solvent has been measured by surface plasmon resonance. 
It is reported that a hydrophilic surface (ho < 0) becomes less hydrophilic upon adding salt, 
whereas no changes have been detected for a hydrophobic surface (h > 0). According 
to Eq. ([2]) the "external" field U for a semi-infinite system (L — > 00) acts like an addi- 
tional, hydrophobic surface field 5h = — J °° zU(z) exp(—z/^) = Aja K^/ (2(1 + k£)) > if 
A7, ctq < 0. A hydrophilic surface becomes less hydrophilic by adding salt or for T — y T c , 
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whereas a hydrophobic surface is influenced only weakly as | tx 1 is small. 

Conclusion. We have demonstrated that even though electrolytes do not alter the uni- 
versal critical behavior of polar solvents close to their critical point, the ion-solvent coupling 
is relevant further away, provided the direct electrostatic interaction is sufficiently weak. 
The crossover from an electrostatics- to a critical Casimir-dominated behavior is expected 
to occur near that temperature at which the bulk correlation length becomes comparable 
with the Debye screening length. Several experiments with monovalent ions in binary liquid 
mixtures can be consistently interpreted in terms of the present model, according to which 
the influence of the ions on the order parameter can be described by an effective "external" 
field proportional to a coupling parameter which measures the difference between the solu- 
bility contrasts of cations and anions in a binary solvent. The insight gained in the present 
study on the effects of ion-solvent coupling may provide an understanding of other situations 
in which critical and non-critical fields are coupled. 
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FIG. 1: Effective wall-wall interaction potential S3 as a function of the scaled wall separation 

kL and the scaled bulk correlation length ac£ for (a) ( ) and (b) (H — ) boundary conditions; 

Csat = k /(^^b) is the saturation surface charge density. For symmetric ( ) boundary conditions 

S3 is repulsive far away from the critical point due to the direct electrostatic interaction between the 
like-charged walls whereas the Casimir force gives rise to an increasing attraction upon approaching 
the critical point. For antisymmetric (H — ) boundary conditions S3 is repulsive far away from as 
well as close to the critical point. Attraction (oc A7) occurs in an intermediate temperature range 
due to the ion-solvent coupling induced by the difference between the solubility contrasts of cations 
and anions in the binary solvent. 
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